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We introduce a new class of models for polymer collapse, given by random walks on regular 
lattices which are weighted according to multiple site visits. A Boltzmann weight uii is assigned to 
each {I + l)-fold visited lattice site, and self- avoidance is incorporated by restricting to a maximal 
number K of visits to any site via setting = for I > K. In this paper we study this model on 
the square and simple cubic lattices for the case K — 3. Moreover, we consider a variant of this 
model, in which we forbid immediate self-reversal of the random walk. We perform simulations for 
random walks up to n = 1024 steps using FlatPERM, a flat histogram stochastic growth algorithm. 
Unexpectedly, we find evidence that the existence of a collapse transition depends sensitively on the 
details of the model. 



I. INTRODUCTION 

The transition of a flexible macromolecular chain from 
a random-coil conformation to a globular compact form, 
called coil-globule transition, has been a subject of exten- 
sive theoretical and experimental studies 1]. Generally, 
polymers in a good solvent are modelled by random walks 
with short-range repulsion (excluded volume). Polymers 
undergoing a coil-globule transition are then modelled by 
adding an additional short-range attraction. The canoni- 
cal lattice model 2, 3] for this transition is given by inter- 
acting self-avoiding walks (IS AW), in which self-avoiding 
random walks on a lattice are weighted according to the 
number of nearest-neighbour contacts. 

From the point of view of continuum models, the draw- 
back of ISAW is that it contains two different kind of in- 
teractions (on-site and nearest-neighbour). In this paper, 
we introduce a different class of lattice models for poly- 
mer collapse, which has only on-site interactions. This is 
in spirit similar to the Domb- Joyce model in which 
a random walk is weighted according to the number of 
multiple visits of lattice sites. 

It is generally assumed that any reasonable random- 
walk model with excluded volume and short-range at- 
traction should describe the coil-globule transition. Ad- 
ditionally, if the collapsed globule is a liquid-like bubble, 
the transition is expected to be second-order 5.], and if 
the collapsed globule is frozen the transition is expected 
to be first-order at least in three dimensions. 

However, an investigation of our new class of mod- 
els reveals that not only the strength of the coil-globule 
transition, but its very existence depends sensitively on 
details of the model. 
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II. THE CLASS OF MODELS AND THE 
ALGORITHM 

We consider 7i-step random walks ^ — (Co, Ci; ■ • ■ i Cn) 
on a lattice. The number of visits to each site x induces 
a density 0j on the lattice sites x via 

n 

Interpreting the density = 0^ as a field induced by a 
particular random walk configuration C, we denote the 
energy of the field as E{(j)). In the Domb- Joyce model, 
the energy functional is given by 

EDjW^aY,m + bJ2^'(^) ■ (2) 

X X 

The first term in this expression is simply related to the 
length n of the random walk, as 

^0(f)=n+l, (3) 

so that a is related to a chemical potential. For 6 = 
we have a pure random walk, while for 6 < the model 
is weakly self-avoiding. The case 6 > leads to an ex- 
tremely collapsed phase, which is dominated by configu- 
rations occupying few lattice sites with very high density. 
Thus, while this model is capable of modelling the swollen 
polymer regime, further terms in the energy functional 
need to be taken into consideration to model "realistic" 
polymer collapse. 

Generalizing Eq. |21 we write the energy for a given 
configuration ^ as 

£;(0 = i?(</'c) = E/('^(^)) • (4) 

X 

In Eq. El /(i) is simply the quadratic polynomial /(t) = 
at -I- \A} , and any particular choice of ][t) gives an alter- 
native to the Domb- Joyce model. 
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Restricting to a maximal number K of visits to any 
site incorporates self-avoidance. Choosing K = \ gives 
self-avoiding walks, and for K > 1 we obtain a model 
with K —\ parameters. To be precise, we choose / to be 
given by /(O) = /(I) = 0, 

/(2) = ei, /(3) = £2, ... f{K)^eK-i (5) 

and f{t) = oo for t > K. Thus, each Z-fold visited site 
contributes to the energy of a configuration. 
The canonical partition function is given by 

where the sum extends over all random walk configura- 
tions with n steps, respectively n + 1 sites. Writing 

e^{£i,...,eK-i) and m = (rrii, . . . , m/^_i) (7) 

where mi denotes the number of sites which are occupied 
by / -|- 1 monomers, the energy can be written as 

K-l 

E{m) — E^ Eirrii — e- fh . (8) 

i=l 

This enables us to write the partition function Eq. as 
Z„(/3) = - (9) 

m m 

where Cn,rh denotes the density of states, and [3 = 
{Pi, . . . , Pk-i) are generalized temperature parameters, 
given by (3i = —f3ei. In other words, {I + l)-fold visited 
sites carry a Boltzmann weight uji = e^' , with t^o = 1 
and uii — for I > K. 

The density of states is estimated directly by the Flat- 
PERM algorithm (see below for a description). Any av- 
eraged quantity Q over the set of parameters m for a 
given length n is calculated by 

{Q)n 0) = . (10) 

fa 

For our simulations, we restrict Xo K — i.e. we only 
allow two-fold and three-fold visits to any site, so that 
we have two free parameters f3\ and 

We consider two variants of the model which differ in 
the underlying set of random walks used. For the first 
variant, we include all simple random walk configura- 
tions, whereas for the second variant, we only include 
simple random walks without immediate self-reversal. 
For this reason, we call the first variant RA for "reversal 
allowed" , and the second variant RF for "reversal forbid- 
den". Clearly, RF configurations form a subset of RA 
configurations. An example of a configuration of the RA 
model is shown in Fig. ^ for the case of a square lattice. 
We shall consider both models in two dimensions on the 
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FIG. 1: Example of a 12-step walk on the square lattice 
with self-reversal allowed (RA). A filled circle (green) denotes 
the presence of a single monomer, filled squares (blue) - two 
and empty squares (white) - three monomers. The numbers 
denote the sequence of monomers. 



square lattice and in three dimensions on the simple cu- 
bic lattice, so that we have a total of four models, which 
we denote by RA2, RA3, RF2, and RF3. 

We have simulated these four models using the Flat- 
PERM algorithm iTj] . The power of this algorithm is the 
ability to sample the density of states uniformly with 
respect to a chosen parametrisation, so that the whole 
parameter range is accessible from one simulation. 

The natural parameters for this problem are m,\ and 
TO2. The algorithm directly estimates the density of 
states Cn,mi,m2 for all n < Umax and any value of mi 
and m2. From this, we can then calculate all interest- 
ing quantities using Eq. ^| As we need to store the full 
density of states, we only perform simulations up to a 
maximal length of Umax = 256. 

Fixing one of the parameters /3i and f32 reduces the size 
the histogram, and enables us to perform simulations of 
larger systems. Fixing (32, say, the algorithm directly 
estimates a partially summed density of states 

Cn,mi {(32) = C'n, 

In this way, we can simulate lengths up to rimax = 1024 
at specifically chosen parameters f3i or (32- Any aver- 
aged quantity ((5)„ is now calculated by using a suitably 
modified version of relation llOl 



III. RESULTS 

For all four models we find SAW behaviour in the case 
of repulsion (i.e. /3i,/32 < 0). Here, singly visited sites 
dominate, and the polymer is swollen, as is clearly ev- 
ident from the scaling of the mean-squared end-to-end 
distance. 

When /32 -C ^ /3i , doubly visited sites should domi- 
nate, and when f3i <C <C /32, triply visited sites should 
dominate. Our simulations confirm this, as well. 
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FIG. 2: Model RF3 with two diflterent phase transitions. On 
varying P2 at fixed negative /3i , there is one type of transition 
(possibly first-order), and on varying /3i at fixed negative P2, 
there is another. The dot represents the point at which the 
type of transition changes. 



We now turn to the question of phase transitions be- 
tween these regimes. Naively one would expect to find 
coil-globule transitions from the swollen phase to the col- 
lapsed region. Moreover, for f3i,P2 ^ 0, there is competi- 
tion between doubly visited and triply visited sites, along 
with the possibility of a further transition. 

We have investigated this scenario in detail for all four 
models. 



A. RF3 

For random walks with forbidden reversal on the sim- 
ple cubic lattice (RF3), we find clear evidence of two 
different phase transitions, leading to the phase diagram 
sketched in Fig. |21 We cannot precisely locate the point 
where the two phase transition lines meet, however, it is 
likely that this point is located in the first quadrant. 

We have analysed these two phase transitions from 
simulations at /9i = —1.0 and (32 — —1.0, respectively. 
Fig. O shows fluctuations in rrii along (32 — -1-0 and 
fluctuations in m2 along (3i — —1.0. In both cases, there 
is a buildup of fluctuations as the system size increases. 
The transition at fixed (32 = —1.0 is stronger than the 
transition at fixed (3i — —1.0. While the latter transition 
is second-order, the former appears to be first-order. It 
may be the case that the latter transition is of the same 
type as IS AW collapse in three dimensions. The first- 
order character of the former transition is supported by 
the fact that the distribution of m2 near the transition 
shows a weak bimodality, see Fig. 0] An investigation 
of the scaling behaviour of the mean-squared end-to-end 
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FIG. 3; Fluctuations in 7712 at /3i = —1.0 (top) and in mi at 
(32 = -1.0 (bottom) for model RF3. 
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FIG. 4; Distribution of m2 at (^2 = —1.0 near the phase 
transition for model RF3. 



distance supports these conclusions. 

There is no indication of any collapse-collapse transi- 
tion in the first quadrant joining up with the point at 
which the type of the collapse transition changes. 
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FIG. 5: Fluctuations in 7712 at /3i — —1.0 (top) and in mi at 
/32 = —1.0 (bottom) for model RA2, showing convergence to 
smooth thermodynamic functions. 



B. RA2 

We now consider random walks with allowed reversal 
on the square lattice (RA2), since it provides the largest 
contrast with RF3. Surprisingly, for RA2, we do not 
find any indication of a phase transition, but merely a 
smooth crossover. Fig.|31shows fluctuations in mi along 
P2 = —1.0 and fluctuations in 7712 along /?i = — 1.0. In 
both cases, there is a smooth crossover, and no buildup 
of fluctuations as the system size increases. There could, 
of course, still be a weak transition. However, an investi- 



gation of the scaling behaviour of the mean-squared end- 
to-end distance supports the conclusion of no transitions. 
At the three points (/?i,/32) = (-1.0,-1.0), (-1.0,1.0), 
and (1.0,-1.0), we find clear evidence for self-avoiding 
walk scaling behaviour. We conclude that RA2 is in the 
self-avoiding walk universality class for all values of (3i 
and ^2- 

So it would seem that changing the dimension and 
allowing for reversals has removed the phase transition 
altogether. This is unexpected. 

C. RA3/RF2 

Our analysis of the two remaining models shows that 
these in some way interpolate between RF3 and RA2. 
Random walks with allowed reversal on the simple cubic 
lattice (RA3) and random walks with forbidden reversal 
on the square lattice (RF2) show behaviour similar to 
each other. 

For negative values of (3i, we find a transition from a 
swollen to a collapsed phase upon increasing /?2. How- 
ever, for negative values of /92, we cannot decide whether 
there exists a very weak phase transition (the specific 
heat exponent a may be negative) or a simple crossover. 
An analysis of the mean-squared end-to-end distance 
scaling is inconclusive. 



IV. CONCLUSION 

In conclusion, we have introduced and simulated var- 
ious new models of polymer collapse in two and three 
dimensions. We have found evidence that the type and 
very existence of the transition depends crucially on sub- 
tle aspects of the underlying lattice model, in particular 
on whether the random walk contains immediate rever- 
sals or not. There is clearly need for further work to be 
done to understand these intriguing results. If backed up, 
these results will surely challenge the current theoretical 
framework of our understanding of polymer collapse. 
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